Bloch electrons in a Jahn- Teller crystal 
and an orbital-density-wave state due to the Berry phase 



00 
OS 

S3 

a 

(N 



> 

m 
o> 
o 

o 

00 



I 

S3 
O 

o 



13 



Hiroyasu Koizumi 

Faculty of Science, Himeji Institute of Technology, Kanaji, Kamigori, Ako-gun, Hyogo 678-12, Japan 

Takashi Hotta and Yasutami Takada 
Institute for Solid State Physics, University of Tokyo, 7-22-1 Roppongi, Minato-ku, Tokyo 106, Japan 

(February 1, 2008) 

The effect of the Berry phase is included explicitly in the wavefunction describing conduction 
electrons in a crystal composed of periodically arrayed Jahn- Teller centers that have conically inter- 
secting potential energy surfaces. The Berry phase can make a drastic change in the band structure, 
leading generally to the formation of an orbital-density-wave state. We discuss implications of our 
theory and possible relations to the orbital ordering observed in the manganese perovskites. 
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The interplay among charge, spin, and orbital degrees 
of freedom is a central issue in recent study on the per- 
ovskite manganites such as Lai^^Ca^MnOa jjj. In par- 
ticular, the direct observation of the orbital ordering ||] 
makes the importance of the orbital degrees of freedom 
even clearer. 

In this work, we investigate Bloch electrons in a Jahn- 
Teller ( JT) crystal by focusing on the effects of double de- 
generacy of e g orbitals at each MnOg octahedron, and its 
lift due to the Jahn- Teller (JT) distortions. Two groups 
have already reported the effects of the JT distortions on 
the electronic properties by implementing first-principles 
band-structure calculations Although these calcu- 

lations are the state of the art, both groups overlook a 
very important ingredient of the problem, namely, the 
Berry phase or the geometric phase [|| associated with 
lattice distortions in each JT center in a crystal. 

In view of this situation, we have done the following: 
First, by starting with a rather general argument on a 
periodic array of JT centers, we specify a simple model 
representing a linear chain of MnOg octahedra. On the 
basis of this model, we describe how the Bloch's theorem 
is modified in the presence of the Berry phase. Second, 
we show that its presence brings about a drastic change in 
the band structure compared to the result in its absence. 
This calls for revision of conventional band-structure cal- 
culations for materials involving the JT centers. Third, 
we find that the optimum JT distortions can be obtained 
by making the total energy minimum and that the resul- 
tant ground state exhibits an oscillation of the orbital 
density, the orbital density wave (ODW). In particular, 
the ODW is reduced to the orbital-ordering state with an 
appropriate choice of parameters involved in the system. 
Fourth, we envisage that many experimental results for 
the perovskite manganites related to the orbital degrees 
of freedom may be explained in the light of our theory. 

Consider a periodic array of JT centers. Here the word 
"JT center" is used in a slightly broader sense than usual. 



We use it for a complex of atoms whose vibronic states 
arise from the interaction of degenerate electronic states 
and linearly-coupled vibrational modes. (Note that the 
degeneracy does not necessarily come from the point- 
group symmetry.) Such a complex has structural insta- 
bility which causes distortions of the atomic configura- 
tion, leading to the lift in the degeneracy of the electronic 
states at equilibrium. If more than two linearly-coupled 
vibrational modes are present, potential energy surfaces 
exhibit so-called "a conical intersection (CI)" ||. The 
CI has also drawn much attention recently, due to the 
topological phase arising from it. This is an example of 
the Berry phase, which has been known in the context of 
dynamical JT effects for years Jt| . 

We specify the JT center we consider in this paper. 
It is composed of two reference electronic states, ip a and 
ipb, and two vibrational modes, Q a and Qb, that couple 
with the electronic states linearly. In matrix representa- 
tion using {ip a , 4>b} as a basis, the vibronic interaction is 
given by 



aQ a f3Qb 



(1) 



This matrix is diagonalized to yield eigenvalues, E = 
{{a + 1 )Q a ±[{a- 1 ) 2 Ql + Ap 2 Ql] 1 / 2 }/2, which possess a 
conically intersecting degeneracy point at Q a = Qb = 0, 
if a =/= 7. 

The isolated JT center described above is modified in 
a crystal. However, the CI seen in E ± survives against a 
perturbation preserving the time-reversal invariance 
Application of such a perturbation modifies the vibronic 
coupling matrix into 



aQ a 

PQb 



a (3Qb 
b jQ a 



(2) 



where real numbers a, 6, and c represent perturbations 
from crystal environments, but this merely shifts the CI 
position to Q a = (c- a) /(a - 7), Qb = —b/fi M. 
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Now we construct a model Hamiltonian H for the crys- 
tal composed of the JT centers. Using the vibronic in- 
teraction expressed in Eq. ([I]), H is given by 

H = -J2(ti a i a j + ^btbj + t 3 a\b 3 + hb\a 3 + h.c.) 

+ y^ y \(aQaj + a)a]aj + (jQ aj + cjb'jbj 
j 

+ (PQkj + b)(a]b J + b}a j )], (3) 

where the first sum describes the electron transfer effects 
between the JT centers and indicates a nearest- 

neighbor pair. The second sum is the electronic Hamil- 
tonian part for the JT center at site j, where Q a j and Qbj 
are Q a and Qb at site j, respectively, and dj and bj are, 
respectively, annihilation operators for electrons in the 
states ip a and tpb- In the later discussion on manganese 
oxides, Oj and bj are, respectively, identified as the elec- 
tronic states steming from d x 2_ y 2 and d 3z i_ r i orbitals of 
Mn 3+ , which are hybridized with oxygen orbitals. 

In Eq. (||), the spin degree of freedom is neglected for 
simplicity, but this is legitimate for the perovskite man- 
ganites due to the Hund's coupling to the electrons in 
the t 2g orbitals. As for lattice vibrations, we have not in- 
cluded the unperturbed phonon Hamiltonian H p h which 
treats the restoring forces for Q aj - and Qbj as well as the 
kinetic-energy parts. In this paper, we restrict ourselves 
to the situation in which the electron density per site is 
constant at all JT centers. Then H p h is reduced to a 
constant and thus it can be neglected. 

Let us transform the electronic basis from {aj, bj} into 
{cj, dj}, where Cj = (ctj+ibj) /\/2 and dj = (aj— ibj)/y2. 
In this new basis, H is rewritten as 

H = - ^(ictcj + td\dj + sc\dj + s*d\ Cj + h.c.) 

(i,3) 

+ YpM c j + 4*i) + v A d J + v ? d H> ( 4 ) 

where t = (h + t 2 )/2, s = (t x - t 2 )/2 + it 3 , Uj = 
[(a + j)Q aj + a + c]/2, and Vj = [(a-j)Q aj +a-c]/2- 
i({3Qbj +b). Note that the circular motion of Q a j and Qbj 
around the CI point {Q a j,Qbj) = (( c ~ a )/( a ~ 7)> — b//3) 
changes the phase of Vj. This phase change is what re- 
sponsible for the appearance of the Berry phase. Because 
we assume that the electron density per site is constant 
at all JT centers, the absolute value of Vj is independent 
of j, though its phase can vary from site to site. 

For the e g electrons in the MnOg octahedron, the sum 
of a and 7 vanishes by symmetry. This makes the param- 
eters Uj constant and these parameters can be removed 
by shifting the origin of energies. For this reason, we will 
not be concerned with the Uj terms. 

We discuss the modification of the Bloch's theorem 
based on the model in Eq. (H) in one spatial dimension 



with t\ = t 2 and t 3 = for the time being. This sim- 
plification leads us to an analytic expression without loss 
of essential physics. In this choice of parameters, H in 
momentum representation is reduced to 



H = Y2 £ k( c l c k + d\d k ) 



k , q 



(5) 



where e k = -2tcosfc and V q = (\V\/N) £V e^^A 
Here we set the lattice constant as unity and N is the 
total number of sites. 

Although it appears in V q , the phase £j does not play 
a role, if it is independent of site, i.e., £j = £. In fact, the 
one-electron eigen energy for H is given by = e^i \V\ 
in this case and the corresponding eigenfunction (p^(x) 
obeys the Bloch's theorem, characterized by the property 



(6) 



with T the operator to translate the whole system by a 
unit lattice vector. 

Imagine now that £j changes from to 2ir site by site 
in the period of M JT centers. For simplicity, we assume 
that £j = 2nj/M and that M is commensurate with N. 
In this situation, V q is given by V q — \V\S 9 2ir/M an d the 
one-electron eigen energy is determined as 



E± = [e k + e k ± V(£fc-e fe ) 2 +4|^| 2 ]/2, 



(7) 



with ik — £k+2ir/M- The corresponding eigenfunction 
<pi(x) satisfies the modified Bloch's theorem: 



9k ( x ) 



where the matrices P and U are defined as 



P 



with pf, 



e itr/M 

e 
given by 





jtt/M 



)' H 








v -Pk 


pi) 



(8) 



(9) 



2 2 N /(£fc-£fc) 2 +4|Vf 



1/2 



(10) 



In addition to the ordinary phase factor e lk , the Bloch 
function acquires the Berry phase. It takes account of 
the difference in the electron motion around the JT cen- 
ter, clockwise or counterclockwise, indicating that careful 
treatments are needed for the analysis of a system includ- 
ing JT centers. 

The phase factor e m l M in Eq. (^J) assures single- 
valuedness of basis functions which satisfy T M ip^(x) = 
e lMk fi (x) in accordance with the Bloch's theorem in a 
crystal with the lattice constant M; without it, signs of 
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basis functions change after the translation by an odd 
M number of sites. This factor may be considered as 
one arising from "a ficticious magnetic field" that exists 
at each CI [[Tof . On the other hand, P accounts for the 
site-to-site variation of the lattice- twisting phase £j . 

The dispersion curves representing are plotted in 
an extended Brillouin zone for various values of M in 
Fig. [I] in which \V\/t is taken to be 0.5 and the origin of 
momenta is shifted by —tt/M in each panel. This shift 
is equivalent to use K = k + tt/M instead of k, where K 
is the generalized quasi- momentum and tt/M arises from 
the ficticious magnetic field | pT| . 

For M — oo (or M = 1), the dispersion relation 
E% = £fc ± |V| is the one corresponding to the result 
in the conventional band-structure calculation [||J|] . For 
M = 2, we observe the same band gap with the magni- 
tude of 2\V\ as that for M = oo, but the dispersion curve 
itself is completely different. We can even predict an in- 
sulating state at "half-filling" , where the electron density 
per site n is unity. For the case of 2 < M < oo, we also 
find a drastic change in E^ from that for M = oo, but 
the system remains metallic irrespective of n. 

For given n, M and |V|/i, information on the dis- 
persion relation allows us to evaluate the ground-state 
energy Em(V) by filling electrons in the one-electron 
states from the bottom. The relative stability among 
the states with different M's can be determined by com- 
paring these Em(V)'s. A typical result is shown for 
\V\/t = 0.5 in Fig. | in which SE(V,M), defined by 
SE(V,M) = E M (V) - E(V = 0), are given as a func- 
tion of n. Note that at V — 0, the ground-state energy 
is independent of M, Thus E(V = 0) is written without 
the suffix M. 

The result in Fig. || indicates that each M has a char- 
acteristic electron filling, n(M), at which 6E(V, M) takes 
a local minimum. For an arbitrary V , n(M) can be eval- 
uated only numerically, but for \V\/t <C I and small 
M, it turns out that n(M) is well approximated by 
n(M) = 2/M (and n{M) = 2-2/M due to the particle- 
hole symmetry). The comparison of SE(V, M)'s suggests 
that the JT-distorted state with the period M is expected 
to occur spontaneously for n around n(M). Note that 
dE/dn exhibits a jump at n = n(M) only for M = 2 due 
to the insulating nature. For other M's, no jump can be 
seen at n = n(M). 

From the viewpoint of our theory, we now try to under- 
stand the orbital ordering observed in Lai-^Ca^MnOs. 
We concentrate our interest into the case of x = 
or LaMnC-3 which corresponds to the half-filling in our 
model. At n = 1, the above discussion on the ground- 
state energy implies that a periodic JT-distorted state 
with M = 2, or an ODW state composed of bipartite A 
and B sublattices, is stabilized. This agrees qualitatively 
well with the experimental situation. 

As for the orbital pattern, LaMnC>3 exhibits an orbital 
order between the orbitals of d 3x 2_ r 2 and d 3y 2_ r 2 rather 



than d x 2_ y 2 and d 3z 2_ r 2. In order to see whether this 
pattern can be predicted in our theory or not, we take 
a more realistic choice of parameters. For this purpose, 
we set the hopping energies as t\ — 3io/4, ti = to/4, 
and tz = v / 3^o/4 by assuming that both ip a and -0b are 
purely <i-wave-like. Here the energy unit to is the over- 
lap integral between adjacent d 3z 2_ r 2 orbitals situated in 
the z-direction. In this choice, t and s are, respectively, 
given by t = to/2 and s — (io/2)exp(«7r/3). On the other 
hand, we set | V| as to by referring to the result of band- 
structure calculations jj| in which the band gap due to 
the JT distortions, 2\V\, is seen to be about the same as 
the bandwidth of the e g electrons. 

With these parameters, we can determine ^ at A and 
B sublattices, £4 and £g, by optimizing the ground- 
state energy, leading to £4 = 7r/3 and £b = 47r/3 (or 
£4 = 47r/3 and £b = 7r/3). The phase difference, 
\£a — £s I = 7T reflects the Berry-phase effect which is 
fixed by the condition M = 2, while the phase 7r/3 takes 
account of the local phase due to the complex number s. 

In order to see the orbital ordering in real space de- 
scribed by the ground state obtained in this way, we 
evaluate the following quantities: 

P3x2-r 2 (i) = {4> x j<t>xj), P3y2-r*(j) = (4>yj<f>yj) , (H) 

where <f) x j and cj> y j are, respectively, the annihilation op- 
erators for electrons in d 3x 2_ r 2 and d 3y 2_ r 2 orbitals at 
site j. These operators are given by the linear combina- 
tion of dj and bj through <f> X j = cos(27r/3)aj+sin(27r/3)&j 
and 4> y j = cos(47r/3)aj + sin(47r/3)6j. 

In Fig. the calculated results for p 3x 2_ r 2 and p^ y 2_ r 2 
are shown as a function of the site index j. The al- 
ternating pattern between d 3x 2_ r 2 and d 3y 2_ r 2 orbitals 
emerges. Thus we may call the present ODW state 
as "the orbital-ordering state". Note that "the orbital- 
ordering state" is obtained without making any fine tun- 
ing of parameters. In fact, the obtained ODW state is 
not easily destroyed by changing the hopping parame- 
ters around the present values. However, it is destroyed 
completely, if we do not pay due attention to the Berry 
phase, as shown in the lower panel in which we take 
U=Cb = 0. 

We have also obtained following speculations: First, 
the interband transitions between E^ and Ef are radi- 
cally different from M = 00 to M = 2. In the former, 
no spectral band is expected but it is expected in the 
latter. This may be connected with the observed anoma- 
lous optical spectra fL2|| . Second, because the Berry 
phase represents a fictitious magnetic field, application 
of "real" magnetic fields provides a profound effect on 
the ODW. This may be the key to explain the observed 
successive phase transitions in magnetic fields |l3| . These 
facts among others may support a Berry-phase scenario 
to understand the orbital ordering in the manganese per- 
ovskites, although we need to make our model more re- 
alistic to compare with the real situations. 
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In conclusion, we have demonstrated the importance 
of the Berry phase in the calculation of band structures 
for crystals containing JT centers. 
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Mitsubishi Foundation. 
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FIG. 1. Dispersion curves for several values of M with 
\V\/t — 0.5. We shift the origin of momenta by —ir/M in 
each panel. 

FIG. 2. Stabilization energy of the M-period JT distorted 
state as a function of the electron filling n for \V\/t = 0.5. 
We note that 5E M {V) is in proportion to \V\ 2 for \V\/t < 1 
and the proportional coefficient depends on M. 



FIG. 3. Orbital densities, p 3l 2_ r 2 and p 3y 2_ r 2, as a func- 
tion of site j for t = to/2 and s — (£o/2)exp(i7r/3). In 
the upper panel, £a and £b are determined by optimizing 
the ground-state energy, while in the lower panel, we set 
U = £b = 0. 
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